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ABSTRACT 

Recurrent nova CI Aquilae has entered the final decline phase a bit before May of 2001, about 300 
days after the optical maximum, showing the slowest evolution among recurrent novae. Based on the 
optically thick wind mass-loss theory of the thermonuclear runaway model, we have estimated the turn- 
off time of the CI Aql 2000 outburst. It is in late March of 2001 after a luminous supersoft X-ray 
source phase lasts ~ 150 days (from November of 2000 until March of 2001). We have also obtained, by 
fitting our theoretical light curves with the 1917 and 2000 outbursts, the white dwarf (WD) mass to be 
M\vd = 1.2 ±0.05 Mq, the helium enrichment of the envelope is He/H~ 0.5 by number, the mass of the 
hydrogen-rich envelope on the WD at the optical maximum is AM max ~ 8.0 x 1O _6 M0, and the average 
mass accretion rate is M acc ~ 1.0 x 10 _7 M Q yr^ 1 during the quiescent phase between the 1917 and 2000 
outbursts. Using these obtained values, we have consistently reproduced the light curve in quiescence as 
well as of the two outbursts. We have also discussed the possibility whether or not CI Aql will explode 
as a Type la supernova in a future. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: individual (CI Aquilae) — 
stars: winds, outflows — X-rays: stars 



1. INTRODUCTION AND SUMMARY 

The second recorded nova outburst of CI Aquilae was 
discovered on 2000 April 28 UT by Takamizawa (2000), 83 
years after the first recorded outburst in 1917 (Rcinmuth 
1925; Williams 2000). CI Aql now becomes a member of 
the recurrent nova class. About 300 days after the opti- 
cal maximum, Hachisu & Kato (2001a) estimated various 
physical parameters of the CI Aql system from the light 
curve fitting and elucidated its nature. They derived the 
white dwarf (WD) mass to be M WD = 1.2 ± 0.05 M Q , 
the helium enrichment of the white dwarf envelope to be 
He/H^ 0.25 by number, and the mass of the hydrogen-rich 
envelope on the white dwarf at the optical maximum to 
be AM max ~ 5.8 x 1O" 6 M . This envelope mass indicates 
an average mass accretion rate of M acc ~ 0.7 x 1O~ 7 M0 
yr _1 during the quiescent phase between the 1917 and 

2000 outbursts. They finally predicted the turn-off time 
of 2001 August and, therefore, that the luminous supersoft 
X-ray source phase lasts until August of 2001. 

After Hachisu & Kato (2001a) has been published, 
some new observational indications appeared concerning 
the turn-off time of the 2000 outburst. The supersoft 
X-ray fluxes were too weak during June and August of 

2001 (Greiner & DiStefano 2002), which indicates that 
the steady hydrogen shell-burning had already vanished. 
Schaefer (2001b) made optical photometry of CI Aql in 



August of 2001 and concluded that CI Aql was V ~ 15.3 
and almost near its quiescent phase (Mennickent & Honey- 
cutt 1995). These results contradict the theoretical pre- 
diction that, if the hydrogen content of the white dwarf 
envelope is A = 0.5 by weight, the steady hydrogen 
shell-burning lasts until August of 2001 (Hachisu & Kato 
2001a). Schaefer (2001c) found on 12 Harvard College 
Observatory archival photographs the 1941 outburst, the 
brightness of which is B ~ 13.8 on May 2 of 1941, and sug- 
gested that the 1917 discovery could indicate a recurrence 
timcscalc of about 20 yrs, with missed outbursts around 
1960 and 1980. In this paper, therefore, we reanalyze the 
light curves of CI Aql, especially reproducing the orbital 
light curves newly obtained in April/May and August of 
2001 (see also Schaefer 2001a,b). 

Additionally, Honeycutt (2002, private communication) 
has revised Mennickent & Honeycutt's (1995) photometric 
data of CI Aql in quiescence. The revised data show about 
0.62 mag down in the U-light curve so that the quiescent 
level is as low as V ~ 16.1, 0.8 mag fainter than August 
of 2001. This value is well consistent with Szkody's (1994) 
V ~ 16.2 in quiescence (see also Table 2 of Mennickent & 
Honeycutt 1995). In this sense, CI Aql did not reach its 
quiescent level but still showed an activity even in August 
of 2001. Therefore, taking these new things into account, 
we have fully reconsidered the CI Aql model. 

To summarize the main revised results, 
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1. The estimated white dwarf mass is the same as the 
previous result, i.e., M\vd = 1-2 ± 0.05 M . 

2. Helium enrichment of the white dwarf envelope is esti- 
mated to be Hc/H^ 0.5 by number, much larger than the 
previous value. 

3. Optically thick winds stopped in early November of 

2000 and the estimated turn-off time is not in August of 

2001 but earlier than May of 2001, i.e., March or April of 
2001. The luminous supersoft X-ray source phase lasted 
from November of 2000 until March or April of 2001. 

4. The mass of the hydrogen-rich envelope on the white 
dwarf at the optical maximum is more massive, i.e., 
AM max ~ 8.0 x 1O~ 6 M than the previous result. This 
envelope mass indicates an average mass accretion rate of 
M acc - 1.0 x 10" 7 M Q yr" 1 between the 1917 and 2000 
outbursts. 

5. These satisfy the conditions of Type la supernova ex- 
plosion if the white dwarf consists of carbon and oxygen. 

In §2, we describe our basic model of CI Aql in the late 
phase of the 2000 outburst and in quiescence. In §3, we 
have theoretically modeled orbital light curves of CI Aql 
in the late phase of the 2000 outburst. In §4, the 1917 out- 
burst is numerically re-analyzed, and we re-determine the 
mass of the white dwarf. In §5, the 2000 outburst is also 
numerically re-analyzed and the duration of the strong, 
optically thick wind and the supersoft X-ray source phase 
are again estimated. Discussion follows in §6. 

2. BASIC MODEL IN LATE PHASE 

Our binary model is illustrated in Figure 1, which con- 
sists of a main-sequence star (MS) filling its Roche lobe, 
a white dwarf (WD) photosphere, and an accretion disk. 
A circular orbit is assumed. We also assume that the sur- 
faces of the white dwarf, the main-sequence companion, 
and the accretion disk emit photons as a blackbody at 
a local temperature which varies with position. For the 
basic structure of the accretion disk, we assume an axi- 
symmctric structure with the size and thickness of 

i?disk = aR*i, (1) 

and 

h=pRdisk {^td ' (2) 

where i?disk is the outer edge of the accretion disk, Rl 
is the effective radius of the inner critical Roche lobe for 
the white dwarf component, h is the height of the surface 
from the equatorial plane, and w is the distance on the 
equatorial plane from the center of the white dwarf. Here, 
we adopt tu-square law (v = 2) to mimic the effect of 
flaring- up at the rim of the accretion disk (e.g., Schandl, 
Meyer-Hofmeister, & Meyer 1997). These two parameters 
of a and (5 are determined by light curve fittings. 

As the August 2001 orbital light curves show non axi- 
symmetric shape, we introduce an asymmetric configura- 
tion of the accretion disk as done by Schandl et al. (1997). 
They explained an asymmetric feature in the orbital light 
curves of the LMC luminous supersoft X-ray source CAL 
87. The cause is due to the gas stream from the com- 
panion hitting the edge of the accretion disk to make a 
vertical spray when the mass transfer rate is as high as in 
the luminous supersoft X-ray sources. We have also intro- 
duced the same type of asymmetry of the accretion disk 



as Schandl et al.'s, i.e., 
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and Clow is a parameter specifying the degree of asymme- 
try, 4>i is the starting angle of the vertical spray, 02 is the 
angle where the vertical spray reaches its maximum height, 
f a is the starting radius ratio from where the vertical spray 
is prominent, z is the height of the disk surface from the 
equatorial plane for the asymmetric case. Here, we assume 
that the accretion disk is symmetric inside w < 0.8i?disk 
{f a = 0.8), and Clow = 0.25, 0i = (17/16)tt, 2 = (9/8)tt, 
unless otherwise specified. Such an example of the accre- 
tion disk is shown in Figure 1. 

The accretion luminosity of the white dwarf (e.g., Star- 
rfield, Sparks, & Shaviv 1988), 
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is introduced, where Lwd is the total luminosity of the 
white dwarf, L\vd,o is t ne intrinsic luminosity of the white 
dwarf, and i?wD is the radius of the white dwarf (e.g., 
i?WD = 0.0068i? Q for 1.2M Q WD). The accretion lumi- 
nosity is as large as 300 L Q for a mass accretion rate of 
M acc - 1 x 10~ 7 M Q yr" 1 onto a 1.2M© WD, although it 
is negligibly smaller than the luminosity of steady hydro- 
gen shell-burning during the outburst, Lwd ~ 26, 000 — 
54,000 L Q , for a 1.2 M Q WD. 

The viscous luminosity of the accretion disk is also in- 
troduced. Then, the original disk surface temperature is 
given by 

3GM WD Macc 

(8) 
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The temperature of the accretion disk surface is as high 
as Tdi s k,org ~ 8500 K at the radius ru = 1 Rq for 
M WD = 1.2 Mq and M acc = 1 x 10~ 7 M Q yr _1 . 

The accretion disk is heated by radiation from a lumi- 
nous white dwarf. Such an irradiation effect is calculated 
as 

L WD a . 3GM WD M a 
Vbk-, ^cos6> + 



disk ' 
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where ?7dk is the efficiency of the irradiation and cos 9 is 
the inclination angle of the surface against light ray from 
the white dwarf (e.g., Schandl et al. 1997). We assume 
?7dk = 0.5, unless otherwise specified. This means that 
50% of the absorbed radiation energy is reemitted from the 
disk surface while the residual 50% is converted to thermal 
energy that eventually heats the white dwarf. The irra- 
diation luminosity (or temperature) is much higher than 
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the accretion luminosity (or temperature) during the hy- 
drogen shell burning, i.e., Tdisk ~ 35,000 K at the ra- 
dius w = 1 Rq for Mwd = 1-2 M Q , cos 9 = 0.1, and 
L WD = 26,000 L Q . 

The outer rim of the accretion disk is not irradiated by 
the white dwarf photosphere so that the temperature of 
the disk rim is simply assumed to be uniform, Tdi s k,rim, 
which is one of the parameter determined by fitting. 

The surface temperature of the irradiated main- 
sequence (MS) companion is roughly estimated by 

<T T MS ~ ^MS^T C0S 9 + ^MS.org- (10) 

where ?7ms is the efficiency of the irradiation, r is the dis- 
tance from the white dwarf, and cos is the inclination 
angle of the surface. The irradiated temperature is as 
high as T MS ~ 26,000 K for r) MS = 0.5, cos (9 = 0.5, 
M WD = 1.2 M (L W d = 26,000 L©), and r - 4i? . 
Therefore, the effect of the irradiation becomes very im- 
portant when the irradiated hemisphere faces toward the 
Earth. 

The temperature of each surface patch is determined by 
calculating all contributions from every patch of the white 
dwarf surface (see Hachisu & Kato 2001b, for more de- 
tail). We have divided the white dwarf surface into 16 
pieces (A6 = n/16) in the latitudinal angle and into 32 
pieces (Acj) = 2n/32) in the longitudinal angle while the 
surface of the main-sequence companion is split into 32 x 64 
pieces, respectively. The surface of the accretion disk is 
divided into 32 pieces logarithmically evenly in the radial 
direction and into 64 pieces evenly in the azimuthal angle. 
The outer edge of the accretion disk is also divided into 64 
pieces in the azimuthal direction and 4 pieces in the verti- 
cal direction by rectangles. When the photosphere of the 
white dwarf becomes very small, e.g., i?disk/-R P h,WD > 10, 
we attribute the first 16 meshes to the outer region (from 

w = i?disk to w = i?disk/\/T0) to avoid coarse meshes 
in the outer part, and then 16 meshes to the inner re- 
gion (from zu = i?disk/v / 10 to w = R p h), each region of 
which is divided logarithmically evenly. The efficiencies of 
irradiation of the main-sequence companion and of the ac- 
cretion disk are assumed to be t]ms = 0.5 and 7/dk = 0.5, 
respectively. 

3. LATE PHASE ORBITAL LIGHT CURVES 

To detect the end of the steady hydrogen shell-burning, 
we have used photometry of CI Aql in May and August of 
2001 (see Figs. 2-4). All photometry was with the 2.1-m 
telescope at McDonald Observatory with the IGI detector 
through a V-filter. The detector is a standard Tektronix 
1026x1024 CCD chip while the filter is the usual Schott 
glass sandwich as prescribed by Bessell (1990) hence re- 
sulting in a close match to the Johnson V-filter. Standard 
IRAF procedures and programs were used to get differ- 
ential photometry of CI Aql with respect to several cali- 
bration stars from Henden & Honeycutt (1995). A total 
of 118 and 429 magnitudes were obtained in April/May 
and August 2001, with a typical photometric uncertainty 
of 0.02 mag. The photometric data are wrapped in one or- 
bital phase, with a new ephemeris obtained by Matsumoto 
et al. (2001), i.e., 

f(HJD) = 2, 451, 701.2086 + 0.61835 x E, (11) 



at eclipse minima. 

The orbital light curve of CI Aql in quiescence were 
obtained by Mennickent & Honeycutt (1995). Recently, 
Honeycutt (2002, private communication) rectified the 
data. Essential points are that their new and old orbital 
light curves are nearly indistinguishable in shape but the 
new data are about 0.6 mag fainter than the old ones. We 
have shown these two light curves in Figures 2-4. 

To fit our theoretical light curves with the orbital pho- 
tometry, we have calculated V light curves by changing 
the parameters of a = 0.5 — 2.0 by 0.1 step, (3 = 0.05 — 
0.50 by 0.05 step, T M s.or g = 5000—10000 K by 100 K step, 
T dis k,rim = 4000—8000 K by 100 K step, L W d,o = 0— 
26,000 L Q by 100 L© step, and i = 70—90° by 1° step 
and seek the best-fit model by the visual comparison of 
the model with the data. After obtaining the best fit 
model for symmetric configurations, we further seek the 
best fit model for August 2000 data by changing the asym- 
metric configuration parameter, Clow, from 1.0 to 0.05 by 
0.05 step. The best fit asymmetric models are also shown 
in Figures 2-4 only for August 2000 data. Dashed lines 
indicate the symmetric case while solid lines denote the 
asymmetric case. The adopted parameters are tabulated 
in Table 1. Enlargement around the eclipse is also shown 
in Figure 5 only for the 1.5M Q MS companion. 

We have calculated the absolute Mb, My, and M/ mag- 
nitudes of our model (see Hachisu & Kato 2001b, in more 
detail). The direct fit with the observational apparent 
magnitude, my, gives us a distance modulus, (m — M)v, 
as shown in Figures 2-4. Since each direct fit determine a 
distance modulus, we have several different distance mod- 
uli for each light curve. We have to choose the most ap- 
propriate distance modulus among them. This value is 
iteratively sought in such a way that one distance mod- 
ulus is consistent all among the light curves of the late 
phase (May/April of 2001), the quiescent phase, the 1917 
outburst, and the 2000 outburst as well as the color excess 
of E(B — V), all of which will be discussed later in more 
detail. 

The light curves are calculated for three companion 
masses of Mms = 1-0 M , 1.5 M and 2.0 M as shown 
in Figures 2-4. Since we have obtained similar light curves 
for these three cases, we cannot determine the compan- 
ion mass without additional information. One may dis- 
tinguish them from the original temperature of the com- 
panion in quiescence; it is too high in the case of 1.0 M 
(7, 500-7, 700 K vs. T ZAM s = 5, 600 K) and too low in the 
case of 2.0 M (6,900 - 6, 700 K vs. T ZA ms = 9, 100 K) 
compared with their zero-age main-sequence (ZAMS) sur- 
face temperature (Bressan et al. 1993). Therefore, we 
adopt Mms = 1-5 M in the following discussion. 

The mass transfer is as high as ~ 1 x 10~ 7 M Q yr _1 , 
and this indicates a thermally unstable mass transfer (van 
den Heuvel et al. 1992), which requires Mms/Mwd ^ 1-1 
(Webbink 1985). Moreover, of our best-fit model, the 
original surface temperature of the main-sequence com- 
panion is roughly consistent with the effective tempera- 
ture of a <~ 1.5 Mq zero-age main-sequence star (TMs.org = 
7,300 - 7, 100 K vs. T ZAM s = 7200 K). Since the radius 
of the 1.5 Mq ZAMS star is about 1.45 i? , this 1.5 M 
main-sequence companion has expanded to fill the Roche 
lobe of <~ 1.7 Rq after the central hydrogen has decreased 
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from X = 0.7 at the zero-age to X — 0.5 at the age of 
~ 1 G yr (Bressan et al. 1993). If the companion has al- 
ready lost significant mass to the white dwarf component, 
its original zero-age mass might be more massive than 1.5 
M . 

The edge of the accretion disk is not fully occulted by 
the companion star at eclipse minima as shown in Figure 
6. This partial occultation implies a shift of the phase po- 
sition at eclipse minima. It may advance when the asym- 
metry of the accretion disk is relatively large and the inner 
side of the spray is brightly illuminated by the white dwarf. 
However, our best fit model in Figure 5 shows no signifi- 
cant advance in the orbital phase at eclipse minima mainly 
because the asymmetry of (j ow = 0.7 is not so large. 

We are able to estimate the turn-off time of steady hy- 
drogen shell-burning on the white dwarf. The orbital light 
curve fitting for May 2001 suggests us that the intrinsic lu- 
minosity of the white dwarf has already decreased to about 
^wd,o ~ 11,000 Lq from 26,000 Lq for the lower limit 
of steady hydrogen shell-burning. The orbital light curve 
fitting for August of 2001 indicates that the intrinsic lu- 
minosity of the white dwarf has much decreased to about 
^wd,o ~ 2, 200 Lq. Thus, we have estimated the turn-off 
time a bit before May of 2001, that is, March or April of 
2001. After the turn-off time, the white dwarf began to 
cool down. It should be noted, however, that hydrogen 
shell burning does not vanish suddenly, but gradually ex- 
tinguishes. Then, our theoretical definition of the turn-off 
time is when the nuclear luminosity decreases below the 
diffusive luminosity of the white dwarf. In other words, 
the diffusive luminosity of the white dwarf is almost bal- 
anced with the nuclear luminosity before the turn-off time. 
Nuclear burning still partly contributes to the luminosity 
of the white dwarf even after the turn-off time. A large 
part of the luminosity in the cooling phase comes from 
the thermal energy in the hydrogen-rich envelope and the 
helium shell (ash of hydrogen burning) of the white dwarf. 

We have calculated the color index to estimate the 
color excess of CI Aql. A rather blue color index of 
(B — V) c w 0.0 is obtained for the case of a 1.5 M© com- 
panion as shown in Figure 7. This suggests a large color 
excess of E(B -V) = (B - V) a - (B - V) c w 1.0 with the 
observed color of (B — V) a <~ 1.0 (see Table 2 of Mennick- 
ent & Honeycutt 1995), where (B — V) c is the calculated 
color and (B — V) is the observed color. This color excess 
is roughly consistent with E(B — V) = 0.85±0.3 estimated 
by Kiss et al. (2001). Therefore, we adopt E(B-V) = 1.0 
in this paper. 

4. THE 1917 OUTBURST 

The light curve of the 1917 outburst has also been re- 
analyzed with the newly obtained system parameters. The 
Tycho B magnitude of the 1917 outburst has been re- 
ported by Williams (2000), which is shown in Figure 8 
together with our calculated B light curves. We have ob- 
tained the apparent distance modulus of (m — M)b = 
15.06 based on the color excess of E(B — V) = 1.0 and the 
^-magnitude distance modulus estimation of (m — M)y = 
14.06. The numerical method for obtaining light curves 
is the same as that in our previous one (Hachisu & Kato 
2001a) and it is fully described in Hachisu & Kato (2001b). 

Assuming a binary system consisting of a white dwarf 



with the mass of Mwd = 1-0, 1.1, 1.15, 1.2, 1.3, and 
1.377 M and a 1.5 M lobe-filling, main-sequence com- 
panion, we have calculated the developing of the outburst 
and obtained B light curves for the hydrogen content of 
X = 0.7. In Figure 8, we do not include the effect of the 
accretion disk or the irradiation effect of the companion. 
These irradiation effects do not contribute much to the 
B light until the B-magnitude decreases to B <~ 12 mag, 
since the white dwarf photosphere has expanded larger 
than the binary orbit (Kato & Hachisu 1994). It is clear 
that there are no mid plateau phases around V ~ 13.5 
for the case of no accretion disk and no irradiation of the 
companion. 

The white dwarf mass is again determined to be Mwd = 
1.2 ± 0.05 M Q from the fitting as seen in Figure 8. The 
various physical parameters of the outbursts are summa- 
rized in Table 2. Here, Z is the content of heavy elements 
by weight, AM max is the hydrogen-rich envelope mass at 
the optical maximum, ?7 w i n d is the ratio of the mass lost by 
winds, AM wind , to AM max , i. = AM wind /AM max , 

t2 is the time taken to drop 2 mag from maximum, tz is 
the time to drop 3 mag from maximum, i w ; n d is the dura- 
tion of optically thick wind phase, and in is the duration 
of steady hydrogen shell-burning. 

We are not able to determine the hydrogen content only 
from the fitting with the 1917 light curve, because much 
lower contents of hydrogen such as X = 0.5 and X = 0.35 
by weight give a similar decline rate or t^) of the early 
phase light curve as listed in Table 2. 

5. THE 2000 OUTBURST 

CI Aql erupted in 2000 April and was densely observed 
in various optical bands (e.g., Kiss et al. 2001; Matsumoto 
et al. 2001). We have determined the hydrogen content 
of the white dwarf envelope from the turn-off time of the 
2000 outburst. 

We summarize the global feature of the optical light 
curve: its optical maximum (V <~ 9 mag) was reached 
on 2000 May 5 (HJD 2451669.5, Kiss et al. 2001). As 
shown in Figure 9, the visual brightness quickly decreased 
to 13.5 mag in about 60 days. Then, it stays at V ~ 14 
mag, i.e., a plateau phase. This mid plateau phase is very 
similar to that of U Sco and can also be explained by the 
irradiation of the accretion disk (Hachisu et al. 2000a, b; 
Thoroughgood et al. 2001). Therefore, we have repro- 
duced the light curve in the plateau phase by assuming 
model parameters similar to those of U Sco. 

Theoretical light curves in Figure 9 are calculated for a 
pair of a 1.2 M white dwarf and a 1.5 M main-sequence 
companion. In the plateau phase, the light curve is deter- 
mined mainly by the irradiations of the accretion disk and 
the main-sequence companion, because the white dwarf 
photosphere becomes much smaller than the binary size. 
The unheated surface temperatures of the companion are 
fixed to be T M s,or g = 7, 300 K throughout the 2000 out- 
burst and Tdi s k,rim = 6, 600 K at the disk rim, which are 
the same as the model parameters in the late phase of the 
2000 outburst. 

The luminosity of the accretion disk depends on both 
the thickness (3 and the size a. Here, we change a by 0.1 
step and (3 by 0.05 step and seek the best fit model for 
the wind phase, the steady hydrogen shell-burning phase, 
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and the cooling phase. Finally we have adopted a set of 
a = 3.0 and — 0.05 during the wind phase to repro- 
duce the orbital modulation observed by Matsumoto et 
al. (2001), a = 0.8 and (3 = 0.30 in the steady hydro- 
gen shell-burning phase, and a = 0.8 and (3 = 0.20 in the 
cooling phase. Here, we adopt symmetric configurations 
of the accretion disk throughout the 2000 outburst. 

The decline of the early phase (t <~ — 70 days) hardly 
depends on the hydrogen content, X, of the white dwarf 
envelope, but the ends of the wind phase and of the hydro- 
gen shell-burning phase depend sensitively on the hydro- 
gen content as tabulated in Table 2. Therefore, we have 
finally determined the hydrogen content of X = 0.35 and 
shown calculated light curves in Figure 9 for three different 
bands of B, V, and I c . 

Here, we have used the apparent distance modulus of 
(m — M)v = 14.06 for the V band. Adopting the color 
excess of E(B — V) = 1.0, we obtain the apparent dis- 
tance moduli of (m — M)b = 15.06 and (m — M)j = 12.46 
for B and I c bands, respectively. Here, we adopted the 
relation of E(I - V)/E(B - V) = -1.60 (e.g., Rieke & 
Lebofsky 1985). Both the B and I c light curves are in 
good agreement with the observational points as shown in 
Figure 9. This indicates the consistency of our estimated 
values. Various physical values of the 2000 outburst are 
summarized in Table 3. 

Matsumoto et al. (2001) reported a sharp <~ 1 mag 
drop of i? c -magnitudc on 23 November 2000 and Kiy- 
ota (2001, VSNET archives, http://vsnet.kusastro.kyoto- 
u. ac.jp/vsnet/) also observed a ~ 1.5 mag drop of I c - 
magnitude around the same day as shown in Figure 9. 
If we attribute this drop to the end of wind phase, the 
hydrogen content of X = 0.35 is consistent with the drop. 

We have re-estimated the luminous supersoft X-ray 
source phase. No luminous supersoft X-rays are expected 
in the early decline phase because the photospheric radius 
of the white dwarf is rather large and the photospheric 
temperature is relatively low. After the recurrent nova en- 
ters the plateau phase, the photospheric temperature of 
the white dwarf increases enough to emit supersoft X-rays 
(Kato 1999). During the massive wind phase, however, 
we do not expect supersoft X-rays because they are self- 
absorbed by the wind itself. U Sco was observed as a 
luminous supersoft X-ray source in the mid plateau phase 
of the 1999 outburst just after the massive wind stopped 
(Kahabka et al. 1999; Hachisu et al. 2000a). Our analysis 
suggests that the massive wind stopped before 23 Novem- 
ber 2000. This means that a luminous supersoft X-ray 
phase started from November of 2000 and continued until 
March of 2001 as shown in Figure 9. The rather low X-ray 
fluxes in June and August of 2001 (Greiner & DiStefano 
2002) are also very consistent with our analysis. 

6. DISCUSSION 

The apparent distance modulus is obtained to be (m — 
M) v = 14.06 by fitting for the late phase of the 2000 out- 
burst. The color excess is also estimated to be E{B — V) = 
1.0 from the difference between the observed color and the 
calculated color, and the absorption of Ay = 3.1 from a 
relation of A v = 3.1 E(B - V) (e.g., Rieke & Lebofsky 
1985). Therefore, the distance to CI Aql is derived to be 
1.55 kpc. The distance modulus of (to — M)b = 15.06 is 
derived from the distance modulus of (to — M)y = 14.06 



and the color excess of E(B — V) = 1.0, which are also 
consistent with the B light curve of the 1917 outburst. 
This set of (to - M) v = 14.06 and E(B - V) = 1.0 are 
also in very good agreement with the V, B, and I c light 
curves of the 2000 outburst as shown in Figure 9. 

The orbital light curve in August of 2001 shows a sig- 
nificant asymmetry. This indicates that the shape of the 
accretion disk is changing in time at least in the final de- 
cay phase of the 2000 outburst. Assuming the asymmetry 
of Clow = 0.7 or 0.75, we have calculated a V light curve 
and added it to Figures 2-4 (solid line for August of 2001). 
On the other hand, the disk shape is almost symmetric for 
May of 2001 and in quiescence. It should be noted that 
the model light curve of August 2001 does not reproduce 
the observational lack of any secondary eclipses even if we 
introduce asymmetry of the accretion disk shape. It may 
stem from large time variations of the edge of the accretion 
disk. 

The envelope mass at the visual maximum is estimated 
to be AM max = 8.0 x 10~ 6 M Q for the hydrogen content 
of X = 0.35, indicating the average mass transfer rate of 
M acc = 1.0 x 1O~ 7 M yr" 1 between the 1917 and 2000 
outbursts, or M acc = 1.3 x 10~ 7 M Q yr _1 between the 
1941 and 2000 outbursts. This mass accretion rate is a 
bit higher but almost consistent with our CI Aql model in 
quiescence (see Figs. 2-4). 

On the other hand, if the mean recurrence period is 
20 yrs as suggested by Schaefer (2001c), the average 
mass transfer rate increases to M acc <~ 4 x 1O~ 7 M 
yr _1 . This value is high enough to maintain steady hy- 
drogen shell-burning for hydrogen content of X = 0.7 
on a 1.2 Mq WD when the mass transfer is steady. In 
other words, it is too high to be compatible with our ther- 
monuclear runaway model, unless the mass transfer rate 
itself rapidly increased just before the outburst. When 
the hydrogen content of the transferred matter is as low 
as X <~ 0.35, the lower limit to steady hydrogen shell- 
burning becomes 8 x 10~ 7 M Q yr _1 (see e.g. eq.(7) in 
Hachisu & Kato 2001b). Then, the mass accretion rate 
of M acc ~4x 10 _7 M Q yr" 1 can still cause shell- flashes 
on a 1.2 Mq WD, although the orbital modulation of the 
light curve in quiescence is not consistent with such a high 
mass transfer rate as shown in Figures 2-4. 

About 82% of the envelope mass is lost by the wind 
(AM wind = 6.6 x 10" 6 M Q ), while the residual 18% 
(AM Ho = 1.4 x 1O~ 6 M ) is left and added to the helium 
layer. The net mass increasing rate of the white dwarf is 
about M Hc = 1.7 x 10~ 8 M Q yr" 1 between the 1917 and 
2000 outburst or M Hc = 2.4 x 10~ 8 M Q yr" 1 between the 
1941 and 2000 outburst. These satisfy the conditions of 
Type la supernova explosion if the white dwarf consists of 
carbon and oxygen (Nomoto & Kondo 1991). Evolution- 
ary paths to Type la supernovae via recurrent novae have 
been discussed in more detail in Hachisu et al. (1999a, b). 

We thank the VSNET members who observed CI Aql 
and R. K. Honeycutt for his providing us the unpublished 
data of CI Aql in quiescence. This research has been sup- 
ported in part by the Grant-in- Aid for Scientific Research 
(11640226) of the Japan Society for the Promotion of Sci- 
ence. 
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Table 1 

Fitted system parameters of CI Aquilae 11 



parameter 


symbol 


M 2 = 1.0 M 


M 2 = 1.5 M 


M 2 = 2.0 M 


inclination angle 


i 


76° 


74° 


72° 


WD luminosity (May of 2001) 




16,000 L Q 


11,000 L e 


7,000 L Q 


WD luminosity (August of 2001) 


^WD,0 


3,500 L e 


2,200 L Q 


1,600 L Q 


MS temperature (May/August of 2001) 


TMS.org 


7,500 K 


7,300 K 


6,900 K 


MS temperature (quiescence) 


TMS.org 


7,700 K 


7,100 K 


6,700 K 


color excess 


E(B - V) 




1.0 




distance modulus 


(m - M) v 




14.06 




absorption 


Av 




3.1 




distance 


d 




1.55 kpc 





The white dwarf mass is assumed to be Mwd — 1.2M© 
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Table 2 

Theoretical decline rates of light curves for the CI Aql 1917 outburst 



1*WD X s * »Jwind b *2 *3 'wind *H 

(M ) (day) (day) (day) (day) 



1.377 


0.7 


0.95 


3.5 


5 


49 


59 


1.3 


0.7 


0.90 


9 


14 


148 


287 


1.2 


0.35 


0.82 


20 


35 


186 


310 


1.2 


0.5 


0.83 


21 


37 


219 


450 


1.2 


0.7 


0.84 


22 


39 


266 


720 


1.15 


0.7 


0.82 


28 


48 


330 


1020 


1.1 


0.7 


0.82 


34 


61 


429 


1450 


1.0 


0.7 


0.82 


55 


78 


730 


1540 



a Solar mctallicity is assumed, i.e., Z — 0.02. 

b The ratio of the mass lost by winds to the accretion mass, i.e., 
AM wind /AM max . 



Table 3 

Theoretical properties of the CI Aql 2000 outburst 



AfwD 


X Z* 


AM max 


AM wind j) wind 


AA/hc 


*2 


*3 


*wi„d 




(Mq) 




(M S ) 


(M ) 


(M ) 


(day) 


(day) 


(day) 


(day) 


1.2 


0.35 0.02 


8.0 x 10~ 6 


6.6 x 10~ 6 0.82 


1.4 x 10~ 6 


18 


30 


180 


307 



a Solar mctallicity is assumed 
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Fig. 1. — Configuration of our CI Aql model in the late phases of the 2000 outburst. The cool component (right) is a slightly evolved 
MS companion (1.5Mq) filling up its inner critical Roche lobe. The north and south polar areas of the cool component are irradiated by 
the hot component (1.2 Mq white dwarf, left). Here, we assume a non-axisymmetric accretion disk in order to reproduce asymmetric light 
curves as shown in Figure 3. The separation is a = 4.25-Rq; the effective radii of the inner critical Roche lobes are R* = 1.53R0, and 
R2 = R2 = 1.69i?0, for the primary white dwarf and the secondary main-sequence companion, respectively. The radius of the white dwarf 
in the cooling phase is as small as i?wD = O.OO68R0, about 1/200 times the Roche lobe size. The WD radius is exaggerated in this figure so 
that it may be seen easily. 
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14 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1.2 M WD + 1.0 M MS (m-M) v =14.06 i = 76° dM/dt=1 .OE-7 

T MS,org = 750 °- 7700 K T disk,rim = 6600 K 0) « = - 8 Z^ ' 15 L WD,0= 1 6OOOL 

(2) a=0.8 £=0.2 asymmetry=0.75 L WDO =35OOL (3) a=0.7 ,8=0.1 L WD0 =0L Q 



(1) 2001 May 



15 



(2) 2001 August 



16 




17 



-0.5 







0.5 



phase 



Fig. 2. — Calculated V light curves are plotted against the binary phase (binary phase is repeated twice from —1.0 to 1.0) together with 
the observational points of the two late phases (filled circles are in April/May of 2001, and small open circles are in August of 2001), and the 
quiescent phase (large open circles: smoothly averaged by us from Mennickcnt & Honcycutt (1995)) as well as the revised quiescent phase 
data (middle size open circles: Honcycutt 2002, private communication). The visual magnitude in quiescence has been revised by 0.62 mag 
down as seen in the figure. Here we adopted the ephemeris given by Matsumoto ct al. (2001) as in equation (11). The model is a binary 
system of 1.2Mq white dwarf (WD) + I.OMq main-sequence (MS) companion. The other model parameters arc shown in the figure and 
Table 1. The solid line for August of 2001 indicates the case of asymmetric structure of accretion disk while the dashed line denotes the case 
of symmetric configuration of accretion disk. The other solid lines correspond to the symmetric accretion disks. Here "asymmetry=0.75" 
means the value of (j ow = 0.75. 
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1.2 M WD + 1.5 M MS (m-M) v =1 4.06 i=74° dM/dt=0.7E-7 

' T MS,org = 730 °- 7100 K T disk,rim = 660 °- 5500 K ( 1 ) « = - 8 ^ = - 15 L WD,0= 1 1 °°°-< 

- (2) a=0.8 £=0.2 asymmetry=0.7 l_ WDO =22OOL (3) a=0.7 (3=0.10 L WD0 =0L Q - 




I I I I I I I I I I I I I I I I I I I 

-1 -0.5 0.5 1 



phase 

Fig. 3. — Same as those in Fig. 2, but for a binary system of 1.2Mq WD + 1.5Mq MS companion. 
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14 



n 1 1 r 



— i 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

1.2 M WD + 2.0 M MS (m-M) v = 1 4.06 i = 72° dM/dt=0.6E-7 

T MS,org = 690 °- 6700 K T disk,rim = 6000 - 5500 K ( 1 ) « = - 8 ^ = - 2 L WDO =7OOOL 

(2) a=0.8 0=0.2 asymmetry=0.7 l_ WDO =16OOL (3) a=0.7 (3=0.1 L WD O =OL 




17 



_i i i i_ 



_i i i i_ 



_i i i i_ 



_i i i i_ 



-1 



-0.5 





phase 

Fig. 4. — Same as those in Fig.2, but for a binary system of 1.2 Mq WD + 2.0 Mq MS companion. 
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i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 p 

1.2 M Q WD + 1.5 M Q MS (m-M) v =1 4.06 i = 74° dM/dt=0.7E-7 

T MS,org = 7300 K T disk,rim = 6600 K « = - 8 ^= - 20 asymmetry=0.7 

L WD0 =2200L Q August 2001 



_l I I l_ 



_l I I l_ 



_l I I l_ 



_l I I l_ 



-0.2 



-0.1 







0.1 



0.2 



phase 



Fig. 5. — Same as those in Fig. 3, but for an enlargement around eclipse minima. Only the data of August 2001 arc shown together with 
the model light curves, solid line for asymmetric and dashed one for symmetric configurations of the accretion disk. 
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Fig. 6.— 



Same as those in Fig.l, but for the configuration at the eclipse minimum. 
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Fig. 7. — Calculated B — V light curve in quiescence is plotted against the orbital phase for a pair of 1.2 Mq WD and 1.5 Mq 
companion. 
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15 



10 - 



m 



15 - 



~ i 1 1 r 



-i 1 1 r 



~ i 1 1 r- 



-i 1 1 r- 



WD + 1.5 M MS i=74° M WD =1.0, 1.1, 1.2, 1.3, 1.377M 
(m-M) B =15.06 X=0.7 Z=0.02 dM/dt=1 .OE-7 T(MS)=7300 K 



1917 outburst of CI Aql 



-i i i L 




300 



350 



400 450 500 

HJD 2421000+ (day) 



550 



600 



Fig. 8. — Calculated B light curves are plotted against time (HJD 2,421,000+) together with the observational points of the 1917 outburst. 
The WD mass is attached to each light curve. Open circles indicate observational points (taken from Williams 2000). The model consists of 
a bloated white dwarf photosphere with no accretion disk and a nonirradiatcd main-sequence companion. The hydrogen content of the white 
dwarf envelope is assumed to be X = 0.70 for all models. The apparent distance modulus of (m — M)b = 15.06 is assumed for all the white 
dwarf masses. Solid lines indicate the light curves connecting the B light at the binary phase 0.5 (roughly the mean brightness) while dashed 
lines correspond to those connecting the B light at the binary phase 0.0 (mid eclipse). 
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10 - 



15 - 



-i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r- 
1.2 IVU WD + 1.5 Mq MS i=74° X=0.35 Z=0.02 dM/dt=0.7E-7 
(m-M) v =14.06 (m-M) B = 15.06 (m-M),= 1 2.46 log R(ph)= 1.40 
T(MS)=7300 T(DK)=6600 a=3.0_0.8_0.8 0=O.O5_O.3O_O.2O 



2000 



tops 




2001 
May 1 



Nov. 1 



ray source phase 



end of steady H — burning 

i 



o o u 



□ 




_l_ 



1600 1700 1800 1900 2000 

HJD 2450000+ (day) 



2100 



2200 



2300 



Fig. 9. — Calculated V, B, and I c light curves are plotted against time (HJD 2,450,000+) together with the observational points of the CI 
Aql 2000 outburst. Small dots (small filled circles, small open circles) indicate observational V and visual magnitudes, while open squares 
represent observational B magnitudes and large open circles indicate observational I c magnitudes (all taken from VSNET archives). Two 
V-magnitudes in the late-phase observed by Schacfer (2001a,b) are indicated by marks. Each light curve connects the brightness at the 
binary phase 0.4. The apparent distance moduli arc (m — M)y = 14.06, (m — M)b = 15.06 and (m — M)i = 12.46 for the V, B, and I c light 
curves, respectively. 



